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ABSTRACT 

We have observed oscillations in the nearby G2 subgiant star (3 Hyi using high-precision velocity 
observations obtained over more than a week with the HARPS and UCLES spectrographs. The os- 
cillation frequencies show a regular comb structure, as expected for solar-like oscillations, but with 
several I = 1 modes being strongly affected by avoided crossings. The data, combined with those we 
obtained five years earlier, allow us to identify 28 oscillation modes. By scaling the large frequency 
separation from the Sun, we measure the mean density of [3 Hyi to an accuracy of 0.6%. The ampli- 
tudes of the oscillations are about 2.5 times solar and the mode lifetime is 2.3 d. A detailed comparison 
of the mixed I = 1 modes with theoretical models should allow a precise estimate of the age of the 
star. 

Subject headings: stars: individual ((3 Hyi) — stars: oscillations 
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1. INTRODUCTION 

There has been tremendous recent progress in ob- 
serving solar-like oscillations in main-sequence and sub- 
giant stars. In a few short years we have moved 
from ambiguous detection s to firm measurements (see 
iBedding fc Kieldsenl [20061 for a recent summary). Most 
of the results have come from high-precision Doppler 
measurements using spectrographs such as CORALIE, 
HARPS, UCLES and UVES. 

The star (3 Hydri (HR 98, HD 2151, HIP 2021) is a 
bright southern G2 subgiant (V — 2.80) that is an ex- 
cellent target for asteroseismology. Early un successful 
attem pts t o measure oscillations wer e made by iFrandsenl 
(fl98l and lEdmonds & Craml <fl995). Five years ago we 
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reported oscillations in th is star from velocity measure- 
ments made with UCLE S (IBedding et al.ll2001l ) and con- 
firmed with CORALIE (jCarrier et all 1200 1[ ). Those ob- 
servations implied a large frequency separation of 56- 
58 /iHz but did not allow the unambiguous identifica- 
tion of individual modes. Meanw hile, theoretical models 
for (3 Hyi h ave been published bvlFernandes k. Monteirol 
(|2003f) and lDi Mauro et all (|2003h . with both studies in- 
dicating the occurrence of avoided crossings for modes 
with I = 1 (see §2]). This goes a long way toward ex- 
plaining the earlier difficulty in mode identification. 

Here we report new observations of (3 Hyi using 
HARPS and UCLES in a dual-site campaign. We con- 
firm the earlier detection of oscillations and are able to 
identify nearly 30 modes, including some which show the 
clear effect of avoided crossings. We also measure the am- 
plitudes of the oscillations and the mode lifetimes, and 
use the large separation to infer the mean stellar density. 

2. VELOCITY OBSERVATIONS AND POWER SPECTRA 

We observed (3 Hyi in 2005 September from two sites. 
At the European Southern Observatory on La Silla in 
Chile we used the HARPS spectrograph (High Accuracy 
Radial velocity Planet Searcher) with the 3.6-m tele- 
scope 12 . A thorium emission lamp was used to calibrate 
the veloc ities, which were extra cted using the HARPS 
pipeline (|Rupprecht et al.1 12004). At Siding Spring Ob- 
servatory in Australia we used UCLES (University Col- 
lege London Echcllc Spectrograph) with the 3.9-m Anglo- 
Australian Telescope (AAT). An iodine absorption cell 
was used to provide a stable wavelength reference, with 
the same setup that we have previously use d with this 
spectrograph for (3 Hyi (|Bedding et al.ll200l] ) and other 
stars. 

With HARPS we obtained 2766 spectra of (3 Hyi, with 
a dead time between exposures of 31 s and exposure times 
of 40 or 50 s, depending on the conditions. With UCLES 

12 Based on observations collected at the European Southern 
Observatory, La Silla, Chile (ESO Programme 75.D-0760(A)). 
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we obtained 1191 spectra, with a dead time between ex- 
posures of 70s and exposure times in the range 60-180 s, 
such that the median sampling time was 187 s (implying 
a Nyquist frequency of about 2.7 mHz). 

The resulting velocities are shown in Fig. [TJ As can be 
seen, the weather was very good in Chile but only moder- 
ately good in Australia (we were allocated 7 nights with 
HARPS and 12 with UCLES). The seeing was also much 
better in Chile. Most of the scatter in the velocities is 
due to oscillations, but there are also slow variations and 
night-to-night variations in both series that we ascribe to 
instrumental effects. Figure [2] shows close-up views dur- 
ing the four brief periods in which both telescopes were 
observing simultaneously. Note that the velocity offsets 
between the two data sets were adjusted separately in 
each segment, to compensate for differential drifts be- 
tween the two instruments. After this was done, we see 
excellent agreement between the two data sets. 

Figure[3]shows the power spectra of the two time series. 
The signal from stellar oscillations appears as the broad 
excess of power centered at lmHz. As usual, we have 
used the measurement uncertainties, Cj, as weights in cal- 
culating these power spectra (according to iu, = 1/of). 
In the case of HARPS, these uncertainties were provided 
by the data processing pipeline. For UCLES, they were 
estimated from the scatter in the residuals duri ng the 
measurement, as described bv lButler et al.l (|1996[) . 

The peak in the HARPS power spectrum at 3070 ^Hz 
(and another at twice this frequency, not shown here) 
is due to a periodi c error in the guiding system , as 
previously noted b y ICarrier fc Eggenbergeri (|2006l ) and 
iBazot et al.l ([2007f ). Fortunately, this signal is restricted 
to a fairly narrow band of frequencies that lies well above 
the oscillations of j3 Hyi, so it does not compromise the 
data. 

In analysing the data, we have followed basically the 
same method that we used for a Cen A dButler et al 
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2005(1 and v Ind (|Bedding et aLll2006f ). Our initial goal 



was to adjust the weights in order to minimize the noise 
level in the Fourier spectrum (§ 12. 1[) . Having done this, 
we then made further adjustments with the aim of re- 
ducing the sidelobes in the spectral window (§ I2.2|) . 

2.1. Optimizing for Signal-to-Noise 

We have chosen to measure the noise in the ampli- 
tude spectrum, cr amp , in two frequency bands on either 
side of the oscillation signal: 230-420 /iHz and 1800- 
2100 /iHz, as indicated by the dotted lines in Fig. [31 We 
averaged these using the geometric mean (since instru- 
mental noise varies as an inverse power of frequency). 
Using this criterion the power spectra in Fig. [3] have 
noise levels (in amplitude) of 6.1 cms -1 for HARPS and 
11 cms -1 for UCLES. These values imply a noise per 
minute of observing time, before any further optimiza- 
tion, of 1.9ms- 1 for HARPS and 3.7ms- 1 for UCLES. 
The difference is due to a combination of factors, primar- 
ily the observing duty cycle and sky conditions. 

The first step in optimizing the weights was to modify 
some of them to account for a small fraction of bad data 
poi nts, in the same way that we have done for other stars 
(see lButler et al.ll2004l for details). In brief, this involved 
(i) cleaning from the time series all power at low frequen- 
cies (below 200 /^Hz), as well as all power from oscilla- 



tions (500-1400 /iHz); and (ii) searching these residuals 
for points that deviated from zero by more than would 
be expected from their uncertainties. We found that 78 
data points from HARPS (2.8%) and 25 from UCLES 
(2.1%) needed to be significantly down- weighted. 

The next step in reducing the noise involved ensur- 
ing that the uncertainties we are using to calculate the 
weights reflect the actual scatter in the data. By this, 
we mean that the estimates of o~i should be consistent 
with the noise level determined from the amplitude spec- 
trum, which mean s they should satisfy equation (3) of 
iButler et al l (|2004ft : 
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We checked this on a night-by-night basis for both in- 
struments, measuring cr amp in the way described above. 
The results showed that the uncertainties u\ for HARPS 
should be multiplied by a factor that ranged from 2.6 
to 4.3, while for UCLES the factor ranged from 0.57 to 
0.73. The variations in this factor from night to night 
presumably reflect changes in the instrumental stabil- 
ity. Note that the factor for UCLES can be compared 
with the value of 0.87 that we estimated for a Cen A 
([Butler et alJl2004t l. which was however determined from 
measuring the noise at much higher frequencies and for 
the run as a whole. For HARPS, the large discrepancy 
between the uncertainties estimated by the pipeline and 
those required to explain the scatter in the data has 
been noted previously, at levels consistent with our re- 
sults (iBouchv et alJ 120051 : ICarrier fc Eggenbergeri 120061 : 
IBazot et al.ll2007| ). 

With the velocity uncertainties corrected in this way, 
we found that HARPS gave a mean precision per spec- 
trum on /3 Hyi of about 1.4 ms- 1 and UCLES gave 
1.7ms -1 (with a spread during the run of about 20%). 
However, HARPS was able to record spectra at about 
twice the rate of UCLES, thanks to the better at- 
mospheric conditions (which allowed shorter exposure 
times) and the faster CCD readout. 

The power spectrum of the combined time series is 
shown in Fig. [4] and the noise level is 4.2 cms -1 in am- 
plitude. We refer to this as the noise-optimized power 
spectrum. Note that the time series has been high-pass 
filtered in order to account for the varying offsets be- 
tween the two data sets and to prevent power from slow 
variations from leaking into the oscillation signal. 

2.2. Optimizing for sidelobes 

The inset in Fig. [4] shows the spectral window (the 
response to a single pure sinusoid) and we see sidelobes 
at ± one cycle per day (11.6 fiHz) that are moderately 
strong (25% in power). A close-up is shown in Fig. [5^. 
These sidelobes occur despite there being relatively few 
gaps in the observing window because of the higher rate 
of data collection with HARPS as compared to UCLES. 

As for our analysis of a Cen A and B, we have also 
generated a power spectrum in which the weights were 
adjusted on a night-by-night basis in order to minimize 
the sidelobes. This involved giving greater weight to the 
UCLES data (multiplying the weights by a factor that 
was typically about 10) and resulted in a spectral win- 
dow with sidelobes reduced to only 4.1% in power (see 
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Fig. [5Jd) . The trade-ofT is an increase in the noise level, 
which rose to 6.2 cms -1 in amplitude. 

Finally, wc also calculated a power spectrum based on 
the best five nights from both sites (JD — 2453600 = 17.5 
to 22.5; see Fig.[lJ. In this case, the weights given to the 
five nights of UCLES data were increased by a single fac- 
tor of 7.2 to lower the sidelobes, whose resulting height 
was 15.7% in power (see Fig. [5};). The noise in this am- 
plitude spectrum was 6.2 cm s -1 . This power spectrum 
has similar noise to the sidelobe-optimized version and 
higher sidelobes, but it has the important property of 
covering a shorter period of time. We expect the modes 
in Hyi to have lifetimes of a few days, and so the best 
signal-to-noise will be achieved in a time series that is 
not too much longer than this. Once the observing time 
greatly exceeds the mode lifetime, the oscillation modes 
in the power spectrum become resolved (eventually split- 
ting into a cluster of peaks under a Lorentzian envelope) 
and the signal-to-noise of the peaks no longer improves 
with increased observing time. 

3. RE- ANALYSIS OF THE 2000 OBSERVATIONS 

To allow us to detect as many oscillation modes as pos- 
sible in Hyi, our analysis included our 2000 June ob- 
servations. These dat a consist of five nights with UCLES 
(jBedding et al.ll200lD and 14 nights with the CORALIE 
spectrograph in Chile ([Carrier et "al"1l200lD , with the lat- 
ter being affected by bad weather. The two time series 
are shown in Fig. [6] 

The UCLES echelle spectra have been completely re- 
processed, using new packages for the raw reduction and 
velocity extraction, resulting in an improvement in preci- 
sion at high frequencies of 13%. We have not re-analyzed 
the CORALIE spectra an d so these velocitie s are the 
same as already published (|Carrier et al.ll200Th . 

We subjected both 2000 velocity time series to the 
same analysis that is described above. The combined 
noise-optimized power spectrum is shown in Fig. [71 
The noise averaged over our defined frequency ranges is 
5.9 cms -1 in amplitude. The spectral window is shown 
in the inset and also in close-up in Fig. [SJl. The sidelobe- 
optimized version of these data has a noise level of 
8.0 cms -1 and its spectral window is shown in Fig. [5^. 

4. OSCILLATION FREQUENCIES 

Mode frequencies for low-degree p-mode oscillations in 
main-sequence stars are well approximated by a regular 
series of peaks, with frequencies given by the following 
asymptotic relation: 

u n ,i = Av(n + fZ + e) - 1(1 + 1)D . (2) 

Here n (the radial order) and I (the angular degree) are 
integers, Aj/ (the large separation) depends on the sound 
travel time across the whole star, Do is sensitive to the 
sound speed near the core and e is se n sitive to the surface 
layers. See IChristensen-Dalsgaardl ([2004T ) for a recent 
review of the theory of solar-like oscillations. 

A subgiant such as Hyi is expected to show sub- 
stantial deviations from the regular comb-like structure 
described by equation ©. This is because some mode 
frequencies, particularly those with I = 1, are shifted 
by avoided crossings with gravity mod es in the stellar 
core (also called 'mode bumping'; see lAizenman et al.l 
Il977f) . These shifted modes are known as 'mixed modes' 



because they have p-mode character near the surface 
but g-mode character in the deep interior. Theoreti- 
cal models of Hyi indeed predict these mixed m odes 
(|Di Mauro et al.ll2003tlFernandes fc Monteiroll2003l ), and 
we must keep this in mind when attempting to identify 
oscillation modes in the power spectrum. The mixed 
modes are rich in information because they probe the 
stellar core and are very sensitive to age, but they com- 
plicate greatly the task of mode identification. 

We have extracted oscillation frequencies from the time 
series using the standard procedure of iterative sine- wave 
fitting. Each step of the iteration involves finding the 
strongest peak in the power spectrum and subtracting 
the corresponding sinusoid from the time series. At each 
step, the frequencies, amplitudes and phases of all ex- 
tracted peaks were adjusted simultaneously to give the 
best fit. In analysing these frequencies, the first step was 
to determine the large separation. We did this by exam- 
ining the pair-wise differences between frequencies ex- 
tracted from the noise-optimized 2005 power spectrum. 
We included all peaks with S/N > 4 (46 peaks) and 
found the strongest signals at 29.1 and 57.2 /iHz. We 
identify these as corresponding to Aj//2 and Av, respec- 
tively, consist ent with the values published from the 20 00 
observations (|Bedding et al.ll2001t [Carrier et al.ll2001h . 

With the large separation established, we next sought 
to identify modes in the echelle diagram. Figure [5] 
shows only the very strongest peaks in the 2005 data 
(S/N > 7). Different symbols are used to show the three 
weighting schemes and the symbol size is proportional to 
the S/N. We can immediately identify the ridges cor- 
responding to modes with I = and 1 = 2, both of 
which are expected to be unaffected by avoided crossings. 
The positions of these two ridges, marked by the vertical 
solid lines, allow us to calculate the expected positions 
of modes with I = 1 and I = 3, based on the asymptotic 
relation (equation^. These are shown by dashed lines 
in the figure. A few of the peaks in that region of the 
diagram fall close to the I = 1 line, but most do not. Wc 
are clearly seeing several mixed modes in Hyi. 

Deciding which peaks correspond to genuine modes is 
always a subjective process. There is a trade-off between 
the desire to find as many oscillation modes as possible 
and the risk of identifying noise peaks as genuine. There 
is also the danger, especially for the weaker peaks, that 
the sidelobe will be stronger than the real peak because 
of interference with the noise, leading to an aliasing error 
of ±11.6 /iHz. And for an evolved star like Hyi, we have 
the added complication of mixed modes. 

In the case of Hyi, we have the advantage of two inde- 
pendent data sets taken five years apart. The 2005 data 
have extremely low noise and a very good window func- 
tion. The 2000 data have higher noise and more gaps, 
but the sidelobe-optimized weighting still gives a good 
spectral window. In Fig. [5] we show all peaks extracted 
from the time series that have S/N > 4. At this level, 
very few of the peaks should be due to noise. We again 
used different symbols to identify the different weighting 
schemes, as shown in the figure legend. 

The circles in Fig. [9] show the final frequencies, which 
are listed in Table [TJ Given the effects of noise and finite 
mode lifetime, we do not expect perfect agreement be- 
tween different measurements of the same mode. In cases 
where a peak was detected in more than one weighting 
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scheme, we have averaged the frequency measurements. 
In those cases, the S/N in the table refers to the maxi- 
mum among the contributors. 

The 28 frequencies above the line are those for which 
we are confident of the identification. These are also 
shown in Fig. Ill) where we see a pattern of mixed modes 
th at is strikingly similar t o that calculated from models 
by iDi Mauro et all ([2003D . Note that the n values in 
Table [T] were determined on the assumption that e falls 
in the range 1 < e < 2, as it does in the Sun. The 29 
frequencies below the line in the table will include some 
genuine modes, some sidelobes that need to be shifted 
by ±11.6 /iHz, and some noise peaks. We list them here 
for completeness, to allow comparison with oscillation 
models of f3 Hyi. 

The uncertainties in the mode frequencies are shown 
in parentheses in Table [1] These depend on the S/N ra- 
tio of the peak and were calculated from the simulations 
described in §(6j assuming the value of the mode life- 
time derived in that section. Of course, these frequency 
uncertainties are based on the assumption that the cor- 
responding peaks in the power spectrum are caused by 
a genuine oscillation modes and are not noise peaks or 
aliases. 

We have looked for a systematic offset between mode 
frequencies of (3 Hyi from the 2005 and 2000 observa- 
tions. Such an offset could indicate variations during a 
stellar activi ty cycle, as has tent atively been suggested 
for a Cen A ([Fletcher et al.ll2006[ ). We find that the 2005 
frequencies are lower than the 2000 frequencies by an av- 
erage amount of 0.1 ± 0.4 /iHz, which is consistent with 
ze ro. Finally, of the tw o mode identifications suggested 
bv lBedding et al. (2001), we can now identify solution B 
as being closest to the correct one. 

4.1. Frequencies from a Bayesian treatment 

We have recently described an alternative approach to 
extracting frequencies from a time series, which we ap- 
plied to the star v In d, using Bayesian methods rather 
than Fourier analysis (|Bedding et al.ll2006t iBrewer et al.l 
I2007T) . We have applied a similar technique to f3 Hyi, us- 
ing the noise-optimized 2005 time series. This analysis 
differed in one important respect from that on v Ind: we 
did not make any assumptions about the distribution of 
frequencies. 

The results are shown in Fig. I10[ where the diamonds 
are the same as in Fig. |H] and the circles are the most 
probable frequencies found by the Bayesian method. In 
most cases there is excellent agreement between the two 
analyses. There is a peak at 1262.20 /iHz, which lies ex- 
actly on the I = 1 ridge, that only appears in Bayesian 
analysis. In fact, this peak was found with S/N = 3.0 by 
the traditional method and we have shown it as a genuine 
mode in Table [1] and Fig. QT] 

4.2. Frequency separations and inferred stellar 
parameters 

The large separation Av is a function of the mode de- 
gree, I, and also varies with frequency. We have deter- 
mined it for both I = and 2 by fitting to the modes 
identified in Table Q] and shown in Fig. QT] The values 
at 1 mHz are given in Table [2] We also give the small 
separation <5z^o2 between / = and 2, and the implied 
values for the parameters Dq and e. 



A detailed comparison of the oscillation frequencies of 
(3 Hyi with theoretical models is beyond the scope of this 
paper. Here, we restrict ourselves to using the large sep- 
aration to estimate the mean density of the star. To a 
good approximation, the large separation is proportional 
to the square root of the mean stellar density: 



Av Q V p Q ' 

However, it is well-known that models - even for the Sun 
- are not yet good enough to reproduce the observed fre- 
quencies due to improper modelling of the surface layers 
(jChristensen-Dalsgaard et all 1 19881 ). There is a system- 
atic offset between observed and computed frequencies 
that increases with frequency, and hence leads to an in- 
correct prediction for the large separation. 

We have examined this surface effect in detail for mod- 
els of the Sun and (3 Hyi (Kjeldsen & Bedding, in prep.) 
and found that the correction term scales in such a way 
that, provided that Av is measured at the peak of the 
oscillation envelope, equation ([3]) remains an excellent 
approximation. We have therefore used equation Q to 
estimate the mean density of (3 Hyi from our measure- 
ment of Av$. We used a value of 134.81 ± 0.09 /iHz for 
the large separation of radial modes in the Sun, which 
we obtained by fitting i n the range n = 17 -25 to the 
frequencies measured by lLazrek et al.l (|1997h using the 
GOLF instrument on the SOHO spacecraft. The result- 
ing mean density for (3 Hyi has an uncertainty of only 
0.6% and is given in Table [21 

5. OSCILLATION AMPLITUDES 

The amplitudes of individual modes are affected by the 
stochastic nature of the excitation and damping. To mea- 
sure the oscillation amplitude of (3 Hyi in a way that is 
independent of these effects, w e have followed the method 
introduced by IKjeldsen et al.l |2005) . In brief, this in- 
volves the following steps: (i) smoothing the power spec- 
trum heavily to produce a single hump of excess power 
that is insensitive to the fact that the oscillation spec- 
trum has discrete peaks; (ii) converting to power density 
by multiplying by the effective length of the observing 
run (which we calculated from the area under the spectral 
window in power); (iii) fitting and subtracting the back- 
ground noise; and (vi) multiplying by A^/3.0 and tak- 
ing the square root, in order to convert t o amplitude per 
oscilla tion mode. For more details, see Kicldsc n et al.l 
<|2005h . 

The result is shown in Fig. [T2] for both the 2005 and 
2000 observations. The difference between the two runs, 
which were made five years apart, is not surprising when 
we recall that the solar amplitude, when measured in the 
same way, shows similar variations from week to week. 
The peak amplitude per mode in (3 Hyi is 40-50 cm s , 
which occurs at a frequency of ^ max = 1.0 mHz. This 
value of !/ max is consistent with that expect ed from scal- 
ing the acoustic cutoff frequency of the Sun ([Brown et al.l 
1 19911: IKjeldsen fc Beddingill995| ). The observed peak am- 
plitude is about 2.5 times the solar value, when the latter 
is measured using stellar techniques (Kjeldsen et al., in 
prep.). Given the uncertainties, this is consistent with 
th e values of 3.2 expected from the L/M scaling proposed 
bv lKieldsen fc Bedding ( 1995]) and also with 2.2 fro m the 
(L/M) - 7 scaling proposed bv lSamadi et all (|2005f) . 
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6. MODE LIFETIMES 

We can estimate the lifetime of solar-like oscillations 
from the scatter of the observed frequencies about th e 
ridges in the echelle diagram (see iKieldsen et al.l |2005() . 
For (3 Hyi we used only the I = and 2 modes, since 
these are apparently unaffected by avoided crossing. In 
addition, we examined the differences between the 2000 
and 2005 data sets for those modes that were detected 
in both (assuming that any differences due to a stel- 
lar activity cycle are negligible; see §2]). We found the 
mean scatter per mode t o be 1.25 ± 0.15 AtHz. Following 
the method described by IKieldsen et all (|2005[ ), this fre- 
quency scatter was calibrated using simulations having 
a range of mode lifetimes and with the observed S/N, 
window function and weights. We assumed that any ro- 
tational splitting could be neglected, which is certainly 
reasonable given the very low measured rotational veloc- 
ity of the star of 2-5kms~ 1 (lDravins fe Nordlundi ri990: 
iReiners fe Schmittl[200l ISetiawan et al.ll2004 T 

We carried out the simulations for ten different values 
of the mode lifetime and the results are shown in Fig- 
ure [T31 where each point represents the mean of 300 sim- 
ulations. The solid curve represents a fit to these points 
and the two dotted curves on either side reflect the 1- 
<7 variations in the simulations. The horizontal dashed 
line shows the observed scatter in the frequencies, from 
which we can read off a value for the mode lifetime of 
2.18lQ; 60 d. 

As a check, we also used a new and completely differ- 
ent method to estimate the mode lifetime. This involved 
examining the region of excess power in the Fourier spec- 
trum and measuring the ratio between the mean in the 
power spectrum and the square of the mean amplitude. 
This ratio, R = (A 2 )/{A) 2 , measures the 'peakiness' 
of the power spectrum. For purely Gaussian noise, R 
will have a value of A/tt = 1.27 (see equation (A2) of 
IKieldsen fc Bedding) Il995f) , while for a spectrum com- 
posed of many strong narrow peaks it will have a much 
higher value. For example, the observed value of R in 
the noise-optimized spectrum from 2005 is 1.42 (mea- 
sured in the range 700-1300 /iHz). Note that the ratio 
of power-to-squared-amplitude is only a sensitive indi- 
cator of mode lifetime in the case where the noise level 
in the spectrum is very low, which is true for our 2005 
observations of (3 Hyi. 

We have again calibrated the observed values of R by 
comparing with simulations with a range of mode life- 
times. The only additional information that we required, 
beyond repeating the assumption of negligible rotation, 
was the number of oscillation modes in the spectrum. 
This is easily calculated from the large separation, pro- 
vided we assume that there is exactly one mode for each 
pair of I and n. For extremely evolved stars in which 
the spectrum of g-modes is very dense, we would expect 
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more than one mixed mode per or der, but for (3 Hyi this 
does not appear to be the case (jDi Mauro et alJ 120031 : 
iFernandes fc Monteirol [20031 ). It is therefore reasonable 
to assume that, despite the presence of mixed modes, 
there is still on average only one mode per order for each 
value of the degree I. 

We again carried out the simulations for ten different 
values of the mode lifetime. The results for the noise- 
optimized spectrum from 2005 are shown in Figure I14( 
where each point represents the mean of 300 simulations. 
The solid curve represents a fit to these points and the 
two dotted curves on either side reflect the l-cr variations 
in the simulations. The horizontal dashed line shows the 
observed value of R, which allows us to determine the 
mode lifetime. We did the same thing for the other two 
weighting schemes for the 2005 data (sidelobe-optimized 
and best five nights). The mean of the three schemes 
gave a mode lifetime of 2.57^J'g2 d. 

These two estimates are independent in the sense that 
they rely on two different properties of the power spec- 
trum. The first measures the positions of the peaks and 
the second measures their amplitudes. We are therefore 
justified in combining the two estimates for the mode life- 
time, and the final value is listed in Table [2] This value, 
combined with the amplitudes estimated in §[5j should 
allow a comparison with the oretical models of the excita- 
tion o f solar-like oscillations (jSamadi et al.1120071 : iHoudekl 
120011 . 

7. CONCLUSIONS 

We have presented dual-site velocity observations of 
j3 Hyi that confirm the presence of solar-like oscillations. 
These data, combined with those we obtained in 2000, 
allowed us to identify 28 oscillation modes. The large 
frequency separation allowed us to infer a very precise 
value for the mean density of (3 Hyi. We also measured 
the amplitudes and lifetimes of the oscillations. Finally, 
the frequencies show the clear signature of mixed I = 1 
modes which, after comparison with models of this star, 
should allow a sensitive measure of its age. 

We would be happy to make the data presented in this 
paper available on request. 
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Fig. 1. — Time series of velocity measurements of Hyi obtained from 2005 August 31 to September 9. 
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Fig. 2. — Close-ups of the time series in Fig.[T]for the four segments during which both telescopes were observing simultaneously. Black 
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Fig. 3. — Power spectra of /3 Hyi from the two instruments, based on the data as it emerged from the pipeline (using the uncertainties 
as weights). The pairs of dotted lines mark the two regions that were use to measure the noise level. 
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Fig. 12. — Amplitude of oscillations in /3 Hyi from the 2000 and 2005 observations (see §|5). 
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Fig. 13. — Calibration of the mode lifetime from measurements of the frequency scatter. Each point shows the mean of 300 simulations, 
each measuring the frequency scatter for a given mode lifetime. The solid curve is a fit to these points and the dotted curves show the 
1-cr variations in the simulations. The horizontal dashed line shows the observed scatter in the frequencies, from which we can read off the 
mode lifetime of j3 Hyi (see §[B}. 
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Fig. 14. — Calibration of the mode lifetime from measurements of the ratio of power-to-squared-amplitude in the noise-optimized power 
spectrum. Each point shows the mean of 300 simulations, each measuring the frequency scatter for a given mode lifetime. The solid curve 
is a fit to these points and the dotted curves show the 1-<t variations in the simulations. The horizontal dashed line shows the observed 
value of the ratio (see §[BJ. 



Bedding et al. 



TABLE 1 
Oscillation frequencies in f3 Hyi 
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TABLE 2 
Parameters for Hyi 



Parameter Value 

Ai/ at 1 mHz (^Hz) 57.24 ± 0.16 

Au 2 at 1 mHz (pHz) 57.52 ± 0.10 

5v 2 at 1 mHz (^Hz) 5.32 ± 0.45 

D at 1 mHz (/xHz) 0.89 ± 0.07 

e at lmHz 1.55 ± 0.05 

p (gem" 3 ) 0.2538 ±0.0015 

mode lifetime (d) 2.32±gf? 



